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1.0 INTRODUCTION 


All spacecraft in low earth orbit are subject to high speed impacts by meteoroids and orbital 
debris particles. These impacts can damage flight-critical systems, which can in turn lead to 
catastrophic failure of the spacecraft. Therefore, the design of a spacecraft for an earth orbiting 
mission must take into account the possibility of such impacts and their effects on the spacecraft 
structure and on all of its exposed subsystem components. 

In addition to threatening the operation of the spacecraft itself, on-orbit impacts also generate 
a significant amount of damaging ricochet ejecta particles. These high speed particles can destroy 
critical external spacecraft subsystems, which in turn also poses a threat to the spacecraft and its 
inhabitants . Ricochet debris particles also increase the contamination of the orbital environment 
and, as a result, constitute a threat to other missions into that environment. Since the majority of 
on-orbit debris inpacts are expected to occur at oblique angles, the characterization of ricochet 
debris created in an orbital debris particle inpact is an issue that must be addressed. 

This report presents a summary of the work performed towards the development of an 
empirical model that that characterizes the secondary ejecta created by a high speed inpact on a 
typical aerospace structural surface. The empirical model developed provides the following 
information as a function of inpact parameters (speed, angle, projectile diameter) and 
target plate geometry (e.g. thickness, etc): 

• angles de fining the spread of ricochet debris and the trajectory of the ricochet debris 
cloud center-of-mass; 

• average velocity of the ricochet debris cloud material; and, 

• velocity and mass of the largest particle(s) in the ricochet debris cloud. 
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In this report. Chapter 2 presents an overview of the phenomenology associated with 
oblique hypervelocity impacts on dun plates, and compares them with the processes typically 
involved in normal (i.e. non-oblique) impacts. Chapter 2 presents a summary of the analysis 
performed to obtain the spatial djs&foutions of ricochet debris particle impacts. This analysis is 
used to determine ricochet debris cloud spray and trajectory angles in terms of impact parameters 
and target plate geometry k 

The technique for calculating the average velocity ofthe ricochet debris cloud is presented 
in Chapter 4. This method is a based on a model developed previously that characterizes the 
masses, trajectories, and velocities of the debris clouds created in an oblique hi^i-speed impact 
[1]. Tlus model employs the thneMoaseivat^ principles, elementary shock physics theory, and 
fundamental thermodynamic principles to obtain a system of algebraic equations for the various 
debris cloud masses, trajectories, and velocities. This emsdng model is modified by incorporating 
the information presented in Chapter 3 and by reducing its dependence on empirical parameters. 

In Chapter 5, relationships for crater diameter and depth are applied to the deepest craters 
in each ricochet witness plate to "back out" the diameters, masses, and velocities of the ricochet 
debris cloud particles that created these craters. These calculations are performed using a method 
similar to that developed in a previous study of ricochet debris particles created in oblique 
hypervelocity impact [2], The information obtained is then used to develop empirical relationships 
that predict the velocity and mass of the largest ricochet debris cloud particle in terms of impact 
parameters and bumper plate thickness. Results obtained using these relationships are compared 
with those obtained previously and presented in Reference [2]. Conclusions derived from the 
work presented herein, as well as recommendations for future activities in this area, are presented 
and discussed in Chapter 6. 
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2.0 OVERVIEW OF HYPERVELOCITY IMPACT PHENOMENOLOGY 
Consider the normal hypervelocity impact of a projectile on the outer bumper of a multi- 
wall system as shown in Figure 2.1. Upon impact, shock waves are set up in the projectile and 
outer bumper materials. The pressures associated with these shocks typically exceed the strengths 
of the materials by several orders of magnitude. For example, in an 8 km/sec aluminum-on- 
aluminum impact, the ratio of the impact pressure (1 16.5 GPa=1.15 MBar) to the strength of the 
material (3 10 MPa for aluminum 6061-T6) is approx. 375, or roughly 2.5 orders of magnitude. 



Pressure Wall P 

Figure 2.1. Hypervelocity Impact of a Generic Multi-Wall System 


As the shock waves propagate, the projectile and outer bumper materials are heated 
adiabatically and non-isentropically. The release of the shock pressures occurs isentropically 
through the action of rarefaction waves that are generated as the shock waves interact with the 
free surfaces of the projectile and the outer bumper. This process leaves the materials in high 
energy states and can cause either or both to fragment, melt or vaporize, depending on the 
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material properties, geometric parameters, and the velocity of impact. 

The outer bumper of the multi-watt structure protects the pressure wall against perforation 
by cttsmtegratiiig fte impacting particle and by creating one or more diffUse debris clouds. In a 
normal impact, only one debris cloud c ontaining both projectile and bun$er plate fragments is 
evident. It first strikes the inner bumper and then travels towards and eventually impacts the 
pressure watt. However, in an oblique impact, three debris clouds are typically formed. Two of 
than, travel inward towards and eventually strike the pressure wall. 

These two debris clouds typically form two distinct areas of damage on the pressure wall. 
Li one damage zone, craters and holes (if any) are nearly circular, which is characteristic of near- 
normal impact. In the other, the craters («d holes, if any) are oblong, indicating that they are 
formed by oblique impacts. As a result, these two debris clouds are often referred to as the 
“normal” and “in-line” debris clouds, respectively. It has hypothesized that the “normal” debris 
cloud contains mainly bumper plate fragments while the “in-line” debris cloud contains mainly 
projectile fragments [3], 

The third debris cloud, often r efe rre d to as the “ricochet” debris cloud, travels backwards, 
away from the multi-wall system. When the projectile obliquity is 45* or less, only a small quantity 
of very fine ricochet debris particles are formed. There can be, however, extensive damage to the 
pressure wall, typically in the form of one or more jagged or petalled holes. As the trajectory 
obliquity is increased beyond 45°, the amount of ricochet debris produced by the impact increases 
significantly. Impacts at obliquities beyond 60* or 65° produce a tremendous amount of ricochet 
debris and only a small quantity of“penetration” debris. The change in behavior that occurs near 
60* has led Schonberg [4] to postulate the existence of a “critical angle of impact obliquity”. For 
aluminum projectiles impacting aluminum bumpers, Schonberg estimated the value of this critical 


angle to be near 60°-65°. Impacts of projectiles with obliquities less than this critical value would 
result in more damage to the pressure wall than to any exterior spacecraft component, while 
impacts at obliquities greater than this critical value would result in more damage to external 
components than to the spacecraft pressure wall. 
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M RICOCHET DEBRIS CLOUD SPREAD AND TRAJECTORY 
In this Chapter, we present a summary of the analyses performed to develop empirical 
equations that define the in-plane spread and trajectory of the ricochet debris cloud in terms of 
impact parameters, material properties and bumper thickness. Tins analysis is based on empirical 
data from two sources: 1)225 high speed impact tests performed at the NASA/Marshall Space 
Flight Center, and, 2) 39 numerical simulation runs performed using SPH, also provided by the 
NASA/Marshall Space Flight Center. 

Figure 3.1 below shows a typical test set-up. This figure is similar to Figure 2.1, except 
that a “ricochet witness plate” has been added to the diagram. These witness plates were typically 
0.3 cm to 1.3 cm thick, depending on the impact conditions, and were provided in each test to 
capture the ricochet debris particles created by oblique impacts. In Figure 3.1, 0 r and 699 denote 
the trajectory of the center-of-mass of the fragments in the ricochet debris cloud and the angle 
below which lies 99% of the damage to the ricochet witness plate, respectively. Based on its 
definition, 899 is presumed to model the spread of the ricochet debris cloud particles. Post-test 
examination of damaged ricochet witness plates revealed several interesting characteristics about 
oblique hypervelocity impact. 


Ricochet 



Figure 3. 1 . Typical Oblique Hypervelocity Impact Test Set-up with Ricochet Debris Cloud 

For impact tests in which the obliquity angle was 30* or less, there was virtually no 
Hamage to the ricochet witness plate. Under such conditions, only a splash deposition was evident 
on that plate. As obliquity increased to 45°, the damage to the ricochet witness plate became more 
pronounced. Small, shallow craters were now evident on the witness plate, typically less than 2 
mm in diameter and less than 2 mm deep, and fairly evenly distributed along the height of the 
witness plate. With further increases in obliquity, an increasing amount of deep cratering became 
evident on the ricochet witness plates. In fact, if a thin ricochet witness plate (i.e. on the order of 
0.3 cm) were used in a test with an obliquity exceeding 65®, it was not unusual to find that the 
witness plate was perforated along the entire length of the border between it and the outer 
bumper. 

From these observations, it became evident that as impact angle increased, the angle 
defining the trajectory of the ricochet debris cloud center-of-mass decreased dramatically, that is, 
as 0 P increased, 0 r decreased. However, even in the high obliquity tests, there were still a fair 
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number of craters near the top of the ricochet witness {dates, indicating that as 6 p increased, 699 
did not experience any significant changes. Appendix A present* a compilation of the 0» and 0, 
data for the oblique impact tests db ft i id ered in this study. Hie value of 0 , for each test was 
obtained by calculating the vertical location of the center of the ricochet witness {date craters 
using a weighted average tpdmli{iie based on die vertical thsirib^on of the witness plate craters. 
The angle 699 was determined based on the height below which lay 99% of the ricochet crater 
damage, and was found simply by counting craters and noting their vertical locations along the 
witness plate. 

To supplement the empirical data, 39 nun^cal runs were performed using SPH, a 
smooth particle hydrodynamics code developed for modelling hypervelocity impact phenomena. 
The impact parameters governing the numerical simulations were chosen to exceed, in terms of 
projectile diameter and impact velocity, those normally attainable with a light gas gun, hi tins 
manner, the “tests” performed using SPH extended the data provided by light gas gun testing. 
Appendix B presents a compilation of the 899 and 6 r data for the oblique impact tests considered 
in this study. For the SPH runs, the value of 0 r for each run was obtained by estimating the angle 
defining the trajectory of the center-of-mass of the ricochet debris cloud based on several SPH 
output plots. The angle 699 was obtained by estimating the angle below which lay 99% of the 
ricochet debris cloud particles as shown on the SPH output plots. 

Three sets of equations for 0 r and 699 were obtained: 1 ) an equations for each based solely 
on empirical data; 2) an equation for each based solely on SPH data; and, 3) an equation for each 
based on a combined database including both empirical and SPH data. These equations are all in 
the following form: 


tan0 r>99 = A| 


'u vrv. 


v d py 


CJ 


COS D 0. 


(3. la-0 


where C* is the bumper material speed of sound. 

Table 3.1 below presents the values of the regression coefficients A-D and the correlation 
coefficients for equations (3.1a-f). Figures 3.2 and 3.3 present a plot of these equations for a 


0.795 cm diameter projectile impacting at 0. 127 cm thick bumper at a velocity of 6.5 km/s at 
trajectory obliquities ranging from 45° to 75°. Also shown in these figures are test data and 
numerical simulation data for 0 r and 099- 

Table 3.1 Parameter Values and Correlation Coefficients for Equations (3.1) 


Equation 

Quantity 

Database 

A 

wm 

C 

D 

Correlation 
Coefficient (R*) 

3.1a 

0r 


0.4725 

0.4085 

0.2299 

0.6458 

0.629 

3.1b 

099 


0.7052 

0.2272 

0.06828 

0.1404 

0.343 

3.1c 

0r 

SPH 

0.1377 

-0.5421 

0.1028 

1.2255 

0.837 

3. Id 

099 

SPH 

1.6519 

0.2201 

0.1689 

1.4587 

0.964 

3.1e 

0r 

Combined 

0.4206 

0.2651 

0.4345 

0.7988 

0.662 

3.1f 

099 

Combined 

0.7608 

0.1989 

0.1146 

0.3191 

0.429 


As can be seen from Table 3. 1, the SPH-only equations have the highest correlation 
coefficients, indicating that the SPH data is very consistent from run to run. In addition, the 
empirical-only and combined equations for 0 r have reasonable R 3 values, which indicates that 
although there is a fair degree of scatter in the empirical 0 r data, the trends in the data are 
consistent over the range of empirical parameters considered. However, as is apparent from the 
very low R 2 values for the empirical and combined 099 equations, there are some features in the 
099 data that are not accounted for in the regression model selected. Additional discussion of theso 
features follow Figures 3.2 and 3.3 below. 
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It is clear from the plots of all three regression equations in Figure 3.2 (empirical-only, 
SPH-only, and combined) that 0 r decreases monotonically as 0 P increases. This is a statistical 
demonstration of the empirical observation made previously regarding the nature of the damage to 
the ricochet witness plates and its relationship to the trajectory obliquity of the impacting 
projectile. 

However, Figure 3.3 shows a divergence in the trend predicted by the SPH-only 
regression and those predicted by the empirical-only and combined egressions. The SPH data and 
the associated curve clearly show a dependence of ©99 on 0 P , one that is similar to that observed 
for 0 r : as 0 P increases, 0<» decreases. However, the empirical data and the associated curves show 
099 to be relatively insensitive to any variation in 0 P . The implication is that the empirical evidence 
dictates that the majority of the ricochet debris cloud particles will always be contained within the 
sflmft spread angle (25° in this case), regardless of the impact parameters. 

The apparent lack of dependence of ©99 on any impact parameter would also explain the 
low correlation coefficients obtained when regressing the 099 data. A multi-variable regression 
process seeks to find trends in the data. When there are none, such as in the case of a constant 
dependent function value, the process returns a correlation coefficient near zero. The discrepancy 
between the empirically-observed independence and the numerically-observed dependence of 099 
is an issue that needs to be explored in more detail in a subsequent investigation. Perhaps more 
consistent calculation (in the case of the test data) and measurement (in the case of the numerical 
data) processes are needed to ensure a more valid joining of the two data sets. 
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4.0 RICOCHET DEBRIS CLOUD VELOCITY 


4.1 Introductory Comments 

A model » developed that can bo used to calculate the masses, vdochities, amt trajectories 
of the three debris clouds created in an oblique hypervelocity impact in toms of impact 
parameters, material properties, and bumper thickness. This model is based on applying the 
principles of mass, momentum, and energy conservation before and after the oblique impact 
event. Etementify shock physics amt thermodynamic principles are used in the model to determine 
the fraction of the initial projectile impact energy that is lost to shock heating of the projectile and 
bumper materials. The model developed is verified by comparing its predictions with available 
experimental information. 

The model is an improvement of the original model developed by Schonberg and Yang [1] 
for two reasons. First, it contains a more widely-applicable empirical equation for 0 r than the 
previous model. Second, it has a decreased dependence on empirical, or user-controlled, 
parameters by explicitly calculating the fraction of the initial projectile kinetic energy that is 
expended in the shock heating and release of the projectile and bumper materials. 

Figure 4. 1 below shows a schematic of the parameters that characterize the motion of the 
three debris clouds created in an oblique hypervelocity impact. In this figure, M|, M 2 , and M, are 
the masses of the ‘normal”, ‘in-line', and ‘ricochet’ debris clouds. Analogously, the quantities Vi, 
V 2 and V„ and 6 t , d 2 , and 6 r are the axial velocities and trajectories, respectively, of these debris 
clouds. We also later introduce the parameter V, (not shown in Figure 4.1) which is used to 
characterize the (assumed equal) radial expansion velocity of each of these three debris clouds. 


01 

02 


M 2 

Figure 4. 1 . Oblique Hypervelocity Impact of a Flat Plate 
4.2 Oblique Impact Model Development 

Applying conservation of momentum before and after the initial impact of the projectile on 
the bumper plate in the vertical and horizontal directions, we arrive at the following equations: 

M P V P cos0 p = Mi Vi cos0, + Mj Vj cos0 2 - Mr V, sin0 r (4.1) 

M P V P sin 0 P = Mi Vi sin 0, + M 2 V 2 sin 0 2 + M, Vr cos0 r (4.2) 

Assuming that no mass is lost in the initial impact, the mass conservation principle yields 
Mp + M f = M| + M 2 + M r (4-3) 

where Mf is the mass of the material that is punched out in the creation of the eliptical hole in the 
bumper plate. This quantity is calculated by noting that for the trajectory obliquities considered, 
the bumper plate hole is elliptical [5]: 

Mf — Jtp b D min D m „tb (4-4) 

where pv and tb are the bumper mass density and thickness, respectively. 

The quantities Dp*, and Dm,* are the lengths of the minor and major axes of the bumper 
plate hole and were calculated using the following empirical equations [4]: 
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0.701 


(4.5) 




cos'" 1 6, +0.93 

0.667 

exp(O.8150, ) + 1.00 


(4.6) 


where C* is the bumper material speed of sound, d. Is foe projectile diameter, and 0, is in radians. 
We note that these equations were derived from hypervelocity impact tests in which spherical 
aluminum projectiles impacted thin aluminum plates. Hence, while the general methodology 
described herein may be valid for other materials besides aluminum, the use of empirical equations 
based on tests employing aluminum plates renders this specific analysis valid only for spherical 
aluminum projectiles impacting aluminum bumper plates. 

Equations (4. 1-4.3) constitute a system of 3 equations in 9 unknowns which must be 
solved for: 3 debris cloud masses, 3 axial velocities, 3 center-of-mass trajectories. An additional 
unknown exists in the form of the average radial expansion velocity of the debris clouds V* which 
must also be solved for. The solution process is facilitated by utilizing experimental observations 
from high-speed impact tests of aluminum dual-wall structures to determine several of the 
unknowns in equations (4. 1-4.3). The remaining unknowns can then be determined in closed 
form. Once this is accomplished, an additional equation can be introduced to solve for V„. The 
process by which this is done is described in the following sections. 

4.3 Trajectory Angles 

The angles 6 t and 6 2 initially increase as 6, is increased [4]. This continues until a critical 
value of0p is reached beyond which 0 t and decrease with continued increases in 8 P . This kind 
of behavior is very difficult to predict analytically without resorting to an advanced shock physics 
analysis. As a result, the analytical prediction of this behavior is beyond the scope of the present 
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work. The empirical equations used to calculate values of 9i and 82 as functions of the initial 
impact parameters are given below [5]: 


®i _ 


e. 


= 0.471 


r 054 

k, J 




'P J 


- 0.086 / \ - 0-478 

' t„ ' 


V*p/ 


cos' 134 e 


cos 0586 0 


(4.7) 


(4.8) 


The angle 9 r is given by the following empirical equation, which was derived in the 
preceding chapter: 


9, =tan -1 



/ > 

t k 

0.2651 

(v\ 

0.4206 



p 


A J 


l c J 


cos 07988 8 


(4.9) 


By using equations (4.7-4.9), 9|, 82 , and 9 r can be treated as known quantities which 
reduces the number of unknowns in equations (4. 1-4.3) to six. 

4.4 Debris Cloud Masses 

The three unknown debris cloud masses are calculated by systematically distributing the 
mass of the projectile and the mass of the bumper plate material that is punched out by the initial 
impact among the three debris clouds and then invoking the conservation of mass equation, 
equation (4.3). This distribution process is accomplished as follows. 

First, it is noted that as 8 P increases, the amount of material in the normal and in-line 
debris clouds monotonically decreases while that in the ricochet debris cloud steadily increases 
[5]. Furthermore, it has been hypothesized that the material in the normal debris cloud is primarily 
bumper plate material, while the material in the in-line debris cloud is primarily projectile material 
[3]. The obliquity of the initial impact on the bumper plate also mandates that the in-line and 
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ricochet debris clouds contain a portion of the bumper plate material. Based on these 
observations, we postulate the following functional forms of Mi and M 2 : 


Mi ■ Mfcos“0 p 

(4.10) 

M a ■ a 2 (Mr - Mf)cos*0p + M p cos*0p 

(4.11) 


where Mf is the mass of bumper plate material that would be qected in a normal impact at a 
reduced velocity V < V„ i.e. Mf = M/(0p=O°, Vp^V), and a 2 is that fraction of the ejected bumper 
{date material in the in-line debris cloud. These forms satisfy the requirement that the debris cloud 
masses decrease as 6, increases and do not violate the hypotheses regarding the origins of the 
material in the respective debris clouds. The values of the exponent n and the coefficient a 2 are 
adjusted so that the final predictions for the debris cloud spread angles based on this analysis 
procedure compare well with those obtained using empirical predictor equations for debris cloud 
spread angles. 

The reduced velocity V used to calculate the mass of bumper plate material in the 'normal' 
debris cloud is taken to be the normal component of the original impact velocity. Any material in 
excess of that which such a normal impact would produce is allocated to the 'in-line' and ricochet 
debris clouds. Therefore, the reduced velocity V is given by 

V' = tiVpCos0p (4.12) 

where t| is a correction factor that is also adjusted so that the final predictions for debris cloud 
spread angles based on the analysis procedure presented herein compare well with those obtained 
using empirical predictor equations. Substitution of equations (4.10-4.1 1) into equation (4.3) 
results in the following expression for the mass of the ricochet debris cloud: 

M, = (1 - a 2 XMf - Mr) cos a 0 p + (M f + M P X1 - cosT0 p ) (4. 13 
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These calculations and assumptions allow Mi, M2, and Mr to be treated as known 
quantities which reduces the number of unknowns to three. Since one of the equations was used 
in the preceding analysis, we now have a system of two equations in three unknowns (V 1, V2, V r ). 
4.5 Debris Cloud Axial Velocities 

Since the 'normal' debris cloud is assumed to contain only bumper plate material and the 
mass of that material is calculated assuming a normal impact, the method for calculating its 
velocity is based on a procedure currently utilized for calculating debris cloud velocities in normal 
impacts of thin plates. This procedure is summarized in the following paragraph. 

The initial normal impact of a projectile on a thin plate produces a shock wave that 
undergoes reflection at the rear surface of the plate. An elementary shock wave propagation 
analysis indicates that the velocity of the rear surface at the moment of reflection is equal to twice 
the particle velocity of the plate material as the shock wave passes through the plate. For a normal 
impact of an aluminum projectile on an aluminum plate, particle velocity is equal to one-half of 
the impact velocity. Hence, a simple substitution shows that for the particular projectile and 
bumper plate materials under consideration, under normal impact, the velocity of the rear surface 
of the plate is equal to the initial normal impact velocity. Since the reflection of the shock wave 
from the rear surface causes the plate material to fragment and thereby creates the debris cloud, 
the presumption is made that the axial velocity of the debris cloud created by the normal impact is 
equal to the velocity of the rear surface of the plate. 

Since the normal velocity assumed to create the 'normal* debris cloud is given by V 1 , then 
the axial velocity of the 'normal' debris cloud is also given by V* , that is, 

V, = tiVpCOse p ( 414 > 

We are now left with a system of two equations in two unknowns, V2 and V r . This system 
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is solved explicitly with the following results: 


M p V p cos(8 p - 9 r ) - V, 008(0, - e r ) 
M 2 cos(02 0 r ) 


(4.15) 


Mp V p sin8 p - Mi V t sin0| - M 2 V 2 sin9 2 
Mf COS0 r 


(4.16) 


Thus, all of the unknowns in equations (4. 1-4.3) are now determined. The final unknown 
to be determined is V* which is found using the method presorted in the next Section. It is 
necessary to determine this unknown in order to be able to validate this model. 

4.6 Debris Cloud Radial Expansion Velocities 

If we apply the principle of energy conservation before and after the initial impact of the 
projectile on the bumper plate, we have the following symbolic equation: 

K.E. U ini-K.E 4 ubm + K.E^ M i (4.17) 

where the initial kinetic energy is that of the incoming projectile, the kinetic energy of the debris 
clouds is that due to their axial motion and expansion, and the kinetic energy that is lost is due to 
the irreversible thermodynamic processes that result from the initial impact such as material 
heating, light flash, etc. If the energy that is lost is written as some fraction £ of the initial impact 
energy, then writing the kinetic energy of the projectile and the debris clouds in standard form 
yields the following: 

|(1 - 5) Mp Vl = |(M, + M 2 + M r ) Vl + |(M, V? + M 2 Vl + M, V?) (4. 18) 


The term on the left hand side of equation (4. 18) may be regarded as the energy available 
for debris cloud motion and expansion. Once the value of £ is known, the only unknown in 
equation (4. 18) is V„ which can be obtained explicitly as follows: 
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fi; ■ 


/ (!-£) Mp v; - (M, Vf + M 2 V| + M, V?) 
V M, + M 2 + M, 


(4.19) 


The parameter which defines the fraction of the initial impact energy that is lost to 
shock heating, is calculated as follows: 




ECM p+ ES'M ( 


2 ’ ’ 


(4.20) 


where E£j and are the waste heats per unit mass produced by the shock heating and 

release of the projectile and bumper hole-out materials. We note that by neglecting energy losses 
such as those due to light flash, the results obtained herein should be conservative in nature. The 
procedure for calculating these waste heats is discussed in the following sub-section. 

4.6.1 Shock Loading and Release Due to High Speed Impact 

In calculating the shock loading and subsequent release of the projectile and outer bumper 
materials, the shock waves are considered to be initially planar. This simplification allows one- 
dimensional relationships to be used for analyzing the creation and release of shock pressures. In 
this manner, the shock pressures, energies, etc., in the projectile and outer bumper materials are 
calculated using the three 1-D shock-jump conditions, a linear relationship between the shock 
wave velocity and particle velocity in each material, and continuity of pressure and velocity at the 
projectile/outer bumper interface. Specifically, if we consider the 1-D impact of a projectile with 
velocity v 0 on a stationary outer bumper, conservation of mass, momentum, and energy across 
the shock fronts in the projectile and in the outer bumper yields 
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u » u »- u » 

V, V,,. 

Frojcctilg: P* -P„ + "V fL (4.21 k) 

V «P 

E„ p =E #p +^(P Hp +P ap XV op -V Hp ) 

u «~ u * 

V,. V* 

u«u_, 

Outer Bumpg: p» = P* +- 77 -^ (4.22a-c) 

V ot 

E»-E.+|(P I ,+P«XV«-V 1 ,) 

where V*l/p is specific volume, u, and ii, are shock and particle velocity, respectively; Vh, Ph, Eh 
and Vo, Po, E, are the density, pressure and energy sties associated with the shocked and initial 
material states, respectively. In equations (4.21a-c) and (4.22a-c), the subscripts 'p', and t' refer to 
projectile and outer bumper quantities, respectively. In the development of equations (4.21a-c) 
and (4.22a-c), the shock velocity in the projectile is taken relative to a 'stationary 1 projectile. 

The linear shock vdodty-particle velocity relationships for the projectile and outer bumper 
materials are taken to be in the form 

u* “ Go klip (4.23) 

where c,=V(KV«) is the material bulk speed of sound, K=E/3(l-2v) is the adiabatic bulk modulus, 
E and v are Young's modulus and Poisson's ratio, respectively, and k is an empirically-derived 
constant. Equations (4.21a-c, 4.22a-c) are applied to the initial impact on the outer bumper of a 
multi-wall system in the following maimer. Upon impact, pressure equilibrium at the projec- 
tile/outer bumper interface implies that 

Pnp = P» (4.24) 
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while material continuity at the interface implies that 

v 0 = Upp + Up, (4.25) 

Because the outer bumper in a multi-wall system is free from any initial mechanical stress 
(it is merely supported at its four comers a fixed distance away from the inner pressure wall), the 
initial conditions ahead of the projectile and outer bumper shock waves are taken to be zero (with 
the exception, of course, of the initial material densities). Solving equations (4.21-4.25) 
simultan eously yields expressions for projectile and outer bumper particle velocities which can 
then be used to calculate shock velocities, pressures, internal energies, and material densities after 
the passage of a shock wave. For example, using this procedure to solve initially for Up, yields 


u 


p* 


b-VA 

2a 


(4.26) 


where 


a = k_ -k, 


'pojJ 

k Pop > 


/ \ 

b = 2k p v o +c 0p +c 0t 

\"op ) 

A = b 2 — 4a(c 0p v 0 +kpV 2 ) 


(4.27a-c) 


Then it follows that 


Upp = v 0 — Up, 

Us, = Co, + k,Up, (4.28a-c) 

Usp = Cop + kpUpp 

The shocked densities of the projectile and outer bumper materials are found by substituting 
equations (4.26, 4.28a-c) into equations (4.21a) and (4.22a) to yield 
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(4.29a) 



u / V 


U_ ” 

■P PP 


Ph« 



“-"Up. 


(4.29b) 


Finally, equations (4.21b,c) and (4.22b,c) are then used to define the pressure and energy 
in the projectile and outer bumper materials, respectively, associated with the passage of the 
shock waves created by the initial impact. Hus completely defines the shocked states of the 
projectile and outer materials due to the initial impact. 

While the shock loading of a material is an irreversible process that results in an increase 
of the internal energy of the shocked material, die release of a shocked material occurs 


isentropically along an 'isentrope' or 'release adiabat'. The difference between the area under the 
isentrope and the energy of the shocked state is the amount of residual energy that remains in the 
material and can cause the material to melt or even vaporize, hi order to calculate the release of 
the projectile and outer bumper materials from their respective shocked states (each characterized 
by Ph, Eh, and V H ), an appropriate equation-of-state is needed for each material. To keep the 
analysis relatively simple, the Mie-Gruneisen equation-of-state [6] was used in this study. 

The Mie-Gruneisen equation-of-state (EOS) is an accurate thermodynamic description of 
most metals in the solid regime and is relatively easy to use. It has the form 

P = P H + pr(E-EH) (4.30) 

where the time-dependent Gruneisen coefficient T is given for most metals as 


r = r ° p » 

P 

In equation (4.31), 


(4.31) 
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(4.32) 


is the ambient Gruneisen coefficient, where K is the adiabatic bulk modulus, p=3a is the 
volumetric coefficient of thermal expansion, and Cp is specific heat at constant pressure. Invoking 


the Second Law of Thermodynamics 
dE = TdS - PdV 


(4.33) 


along with the isentropic constraint dS-0 for the release process allows us to construct the release 
isentrope in P-V space for a material referenced to the material Hugoniot in P- V space and a 


given initial shocked state defined by P H , V H , Eh. Using the procedure outlined in Reference [6], 


the pressure Pi at a specific position *i* along the isentrope can be shown to be given by 



(4.34) 


where AV is the incremental change in volume used to create the release isentrope, and Phi and 
Em are the pressure and energy along the Hugoniot corresponding to the i-th position in the 
release process. The release process is continued using equation (4.34) until the release isentrope 
so determined crosses the V-axis. 

It should be noted that based on its formulation, the Mie-Gruneisen EOS cannot be 
expected to give accurate results in a highly expanded liquid regime or in a vapor regime. This is 
because as impact energy increases, the assumption that the Gruneisen coefficient is a function of 
density alone is no longer valid. At high impact energies, the Gruneisen coefficient is a function of 
internal energy as well as density. Experience has shown, however, that it does yield fairly ac- 
curate end-state results even when there is a small percentage of molten material present [7]. 
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Once the release process calculations for the projectile and bumper materials have been 
completed, the areas under the respective isentropes are calculated and subtracted from the initial 
shocked energy state to determine the respective waste heats, that is, 

= e h “ A. (4.35a) 

= e h ~ A-kwT (4.35b) 

4.7 Oblique Impact Model Vafficatea 

The validity of the proposed method of solution for the ten unknowns that characterize the 


the attendant assumptions) is assessed by comparing model predictions of debris cloud spread 
angles with the predictions of empirically based equations for debris cloud spread angles. Model 
values for the spread angles of the 'normal* and 'in-line' debris clouds, and respectively, are 


given by: 

♦i=2tan''(^j I - 1, 2 


(4.36) 


The empirical values of debris cloud spread angles are found using the following relationships [5]: 


/ v \ a 90 Y# V 

It ) cos °”‘ 9 


1.906/ \ 0.345 


cos 0 ™ 0 


(4.37a) 


(4.37b) 


Table 4.1 presents the a summary of the impact paramters used in the evaluation of the 
model developed herein. Tables 4.2a-c, 4.3a-c, and 4.4a-c present foe final values of the user- 
controlled parameters 012 , q and n corresponding to the impact conditions in Table 4.1. 
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Table 4.1. Impact Conditions Considered in Model Validation 





Values Considered 


4.0. 5.5. 7.0 


30, 45,60 


0.635. 0.795, 0.953, 1.13, 1.27 


1.3. 1.6. 2.0 


Bumper Thickness, tb (mm 


Table 4.2a. Model Parameters aa, r| and n for 0,=3O°, tb-1.3 mm 



N 

aa 

3.45 

1.00 

2.40 

1.00 

1.50 

1.00 

0.35 

1.00 

3.45 

1.00 

2.45 

1.00 

1.40 

1.00 

0.60 

1.00 

3.40 

0.95 

2.50 

0.93 

1.50 

0.91 

0.90 

0.89 



Table 4.2b. Model Parameters ot 2 , n and n for 0,-30® 4-1.6 mm 



n 


4.0 i 

0.635 

0.85 

4.50 

am 

0.795 

0.95 

3.40 

4.0 

0.953 

1.05 

2.50 

4.0 

1.13 

1.15 

1.60 

5.5 

0.635 

0.75 

4.50 

5.5 

0.795 

0.85 

3.45 

5.5 

0.953 

0.95 

2.60 

5.5 

1.13 

1.05 

1.80 

7.0 

0.635 

0.75 

4.40 

7.0 

0.795 

0.80 

3.50 

CM 

0.953 

0.85 

2.70 

7.0 

1.13 

0.90 

1.90 
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Table 4.3a. Modd Parameters 0 % r], rad n for 9 P =45°, 4 = 1 .3 mm * 


<*a 


4.0 


1.00 

1.85 

BSiSi 

MXM 

0.795 

1.10 

1.35 

1.00 

4.0 

0.953 

1.35 

0.85 

LOO 

WSM 

111 

1.50 


1.00 



0.95 1.95 1.00 


mam 


1.05 . 


1.00 

5.5 

1.10 

0.85 

LOO 

5.5 

1.13 

1.15 

gyiin 


0.635 

0.85 

1.85 

0.95 

IKZB 

0.795 

0.95 

1.35 

0.93 


0.953 

1.05 

0.90 

0.91 

WBM 

1.13 

BWK1 

0.45 

0.89 
































































Table 4.3b. Model Parameters 012 , tj, and n for 0 P =45°, 4=1.6 mm 


m 

mi 

D 

■ 

eta 

mem 

0.635 

1.00 

2.45 

1.00 

MniMt 

0.795 

1.10 

1.90 

1.00 

mm 

0.953 

1.20 

1.40 

1.00 

MEM 

1.13 

1.35 

1.00 

1.00 

1 5.5 

0.635 

0.95 

2.55 

1.00 

5.5 

0.795 

1.05 

1.95 

1.00 

5.5 

0.953 

1.10 

1.45 

1.00 

5.5 

1.13 

1.15 

1.05 

1.00 

mem 

0.635 

0.85 

2.45 

0.95 

mem 

0.795 

0.95 

1.90 

0.93 

1 7.0 

0.953 

1.05 

1.45 

0.91 

mem 

1.13 

1.15 

1.05 

0.89 


Table 4.3c. Model Parameters a*, rj, and n for 0 P =45°, 4=2.0 mm 


n 

KH 

n 

ma 

a 2 


kgCjI 

M 



4.0 

0.635 

1.00 

3.05 

1.00 

4.0 

0.795 

1.10 

2.50 

1.00 

4.0 

0.953 

1.20 

2.00 

1.00 

4.0 

1.13 

1.35 

1.60 

1.00 

5.5 

0.635 

0.95 

3.05 

1.00 

5.5 

0.795 

1.05 

2.45 

1.00 

5.5 

0.953 

1.10 

2:00 

1.00 

5.5 

1.13 

1.15 

1.65 

1.00 

7.0 

0.635 

0.80 

2.85 

0.95 

\mem 

0.795 

0.90 

2.40 

0.93 

7.0 

0.953 

1.00 

1.95 

0.91 

7.0 

1.13 

1.10 

1.65 

0.89 
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Table 4.4a. Model Parameters ct 2 , r|, and n for 0^*60*, 4=1.3 mm 
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Table 4.4c. Model Parameters a 2 , r| and n for 0 P =6O°, 4=2.0 mm 


m 

mm 

n 

n 

a> 

TIThT/II 

■SSI 

Wm 



KOI 

0.635 

1.50 

2.25 

1.00 

KXB 

0.795 

1.50 

1.90 

1.00 

mm 

0.953 

1.55 

1.60 

1.00 

HEEH 

1.13 

1.75 

1.35 

1.00 

5.5 

0.635 

1.30 

2.20 

1.00 

5.5 

0.795 

1.40 

1.90 

1.00 

5.5 

0.953 

1.50 

1.60 

1.00 

5.5 

1.13 

1.60 

1.40 

1.00 

\MBEM 

0.635 

1.15 

2.10 

0.95 

7.0 

0.795 

1.25 

1.80 

0.93 

mm 

0.953 

1.35 

1.60 

0.91 

HEZ3H 

1.13 

1.45 

1.40 

0.89 


F inall y Table 4.5a-c present percent error summaries showing differences between 
prediction and experiment for the various bumper plate thicknesses, impact trajectories, projectile 
diameters, and obliquities considered. For each perforating debris cloud spread angle, the value 
shown is the precent difference between model prediction and empirical equation prediction. As 
ran be seen from Table 4.5a-c, the values of the spread angles that result from the calculations 
described herein are very close to the experimental values. Naturally, the values of the parameters 
a 2 , T] and n have been adjusted to ensure that model predictions and empirical results are closely 
matched. 
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Table 4.5a. Percent Error Summaries for U = 1.3 mm 


1 V. - 4.0 km/s | 

«*p 

(cm) 

i 30 dec 

45 det 

L IP det 1 

rn 

rsz 

mm 

mm 

mm 


0.635 

IMFEff 

■m 


WM\ 


0.64 

0.795 

1 -1L76 

11.01 

WBlEl 

mm 

1 -9.23 

3.48 

0.953 

IKRUR 

22.52 


ffESEi 

mmi 

16.24 

1.13 

\mm 

55.16 

WzSSl 

■7f!T 

1 -18.35 

19.79 

1 V.-5.5km/« 1 


1 30 det 

45 det 

I 60 det 1 

\MSSM 

Mm 

n 

mm 

b 

h 

4a 

0.635 

1.19 


MCSI 

4.25 

4.59 

8.05 

0.795 

mm 

WEMl 

wem 


WESB 

3.41 

0.953 

Mtm 

■rain 

mm 


-7.89 

14.63 

1.13 



SEES 

smt 

BUS 

20.57 

■ 1 

1 «. 

30 det 

45 det 

60 det 1 

KS5B 

4l 

4a 

■m 

h 

mm 

KB 

0.635 

-0.30 

3.63 

8.57 

6.22 

2.09 

3.79 

0.795 

BEE2 

-3.27 

2.11 

4.94 

-3.81 

-1.01 

0.953 

-6.95 

mm 

-2.69 

11.45 

WEB® 

8.05 

1.13 

-26.23 

14.18] 

-4.24 

26.97 

-9.86| 

14.57 


30 


iiii- 


B=fc-H N=Fa4i-y-^ ■ 
























































Table 4.5b. Percent Error Summaries for t* = 1.6 mm 


1 V- * 4.0 km/s I 

mm 

HZ 

leg 

45 dec 

60 dec 1 

■&3I 

mm 

MM 

mm 


n 


0.635 

-1.30 

4.97 

1.83 

2.20 

- 3.59 

-3.721 

0.795 

-9.38 

7.50 

-7.50 

1.26 

4.43 


0.953 

-15.55 

15.73 

-8.33 

16.09 

4.96 

18.931 

1.13 

-17.19 

35.45 

-16.88 

27.99 

-15.71 

St 

00 

I V. - 5.5 km/s 1 

mm 

30 dec 

45 dec 

60 dec 1 

■551 

mm 

MM 

MM 

MM 

MM 

h 

0.635 

1.23 

3.83 

-3.15 

0.62 

-2.16 

1.74 

0.795 

-5.85 

3.79 

-8.65 

3.54 

-4.79 

0.18 

0.953 

-13.02 

7.67 

-6.42 

13.20 

-5.48 

9.39 

1.13 

-18.37 

16.95 

-9.89 

20.85 

-13.71 

12.97 

1 V B * 7.0 km/s 1 

mm 

■m 

45 dec 

60 des 1 

■SSI 

n 


mm 

MM 

MM 


0.635 

-8.26 

9.55 

-2.56 

1.72 

-5.71 

5.12 

0.795 

-9.06 

-1.28 

-4.63 

2.66 

-5.92 

1.94 

0.953 

-11.65 

-3.94 

-8.89 

6.60 

-6.00 

10.71 

1.13 

-12.07 

1.43] 

-13.33 

13.38 

-8.95 

17.75 
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Table 4.5c. Percent Error Summaries for t* = 2.0 mm 


v. - 4.0 tart __ 1 

* 

(««) 

1 30* 

is 

1 45 dec 

L <0dec 1 

ra 

ra 

m 

m 

i 

\wm 

EEa 

MEiVvl 

lags 

IHSEE 

Hfe: 

mm 

mm 



HEE3 

iBEES 

KHi 

mm 

MM 

Bfgi 1 

BE 


WBSR 

id 

mm 

HE 

mwm 

Mmnmni 

BIS 

ra 


18.12 


BE 

l mh i 

* 

L »«fer 

45 1 

is 

l. : 60 dec 1 

n 

n 

♦l : 

L*a_ 

mm 

I 

roi 


9.13 

iP 1^. 

BO 

BEE 



■ESB 

9.67 

flTTi 

BE 

MEM 

1.18 

0.953 

-4.4* 

935 

-10.13 

13.93 

KH 

9.91 

1.13 



BHit M 

13.79 

BEE 

10.00 

f V.- 7.0 tart I 

dr 

-(«■> 

30 dec 

45 1 

!s J 

60 deg I 

B&BB 

il 

n 

mm 

IH 

fc 

EZSXi 

BSK8£i- 


Has 

anm 


7.85 

0.795 

-7.69 

11.44 


mm 

WEM 

8.47 

0.953 

-7.80 

6.33 

ihi 

KEE2! 

-4.13 

0.19 

1.13 

HEOI 

9.03 

-16.631 

5.3 1| 

-9.58 

-3.26 





























5.0 CHARACTERIZING RICOCHET DEBRIS CLOUD PARTICLES 
Damage potential estimates of ricochet debris particles created in an oblique hypervelocity 
impact will contribute significantly to the successful design of an effective protection systems for 
external spacecraft components and will assist in determining the overall survivability probability 
of a spacecraft following such an impact. A simple way of modelling the damage potential of a 
ricochet debris particle is through its size and speed. 

In this Chapter, a technique is presented for developing empirical relationships that predict 
the velocity and mass of the largest ricochet debris cloud particle in terms of impact parameters 
anH bumper plate thickness. This is accomplished by "backing out" the diameters, masses, and 
velocities of the ricochet debris cloud particles from measured craters penetration depths and 
surface diameters on the ricochet witness plates of 139 oblique impact tests performed at the 
NASA/Marshall Space Flight Center. Measured values of crater depth and diameter are used 
together with empirical relationships for these quantities to determine particle diameters and 
velocities. Results obtained using these relationships are compared with those obtained previously 
and presented in Reference [2]. Visual inspection of damaged ricochet witness plates reveal 
several interesting features that address the validity of this method. 

1) The surface openings of ricochet witness plate craters formed by debris impacts were 
very nearly circular, which is indicative of near-normal impact trajectories. This 
observation is confirmed by the analysis performed in the preceding chapter, which 
concluded that most of the ricochet debris particles will be contained within a cone 
having an apex angle of 30o or less, regardless of the original impact angle. 
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2) In the tests where the ricochet witness plates were sufficiently thick, the reverse sides 
of the plates remained smooth and undamaged even though the front sides exhibited 
si&iific«Fit crater damage. In these cases, the post-impact appearance of the ricochet 
witness plate was identical to that of a "thick plate" subjected to the same debris 
particle impact loading. 

Based on these observations, the use of thick plate equations for penetration depth and 
crater diameter due to normal hypervelocity impact is justified provided that the reverse side of 
the ricochet witness plate in which the crater depths are measured is smooth and undamaged (i.e., 
no spall or dunphiig). 

Examination of existing penetration depth equations revealed a strong coupling between 
particle size and velocity effects. That is, the same size crater can be produced by a small particle 
traveling at a high speed or by a larger particle traveling at a slower speed. Therefore, in order to 
have a unique solution for particle size and speed, a second set of equations describing another 
measurable crater quantity was needed. A search of existing literature on cratering phenomena in 
hypervelocity impact suggested crater volume to be such a quantity. Thus, a crater volume 
equation used in conjunction with an equation for penetration depth could be used to solve 
uniquely for particle size and speed. Since it is more facile to measure the surface itiamatw of an 
impact crater than it is to determine its exact volume, the crater volume equations were rewritten 
in terms of surface diameter. The analysis dim proceeded as follows. 

First, penetration depths and surface diameters of the three hugest craters on ricochet 
witness plates with undamaged rear surfaces were measured ((dates with through-holes or only 
splash damage were not considered). Second, crater volumes were calculated for each measured 
crater. The crater with the largest volume, deemed the most damage as a result, was identified and 
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retained for future analysis. By considering only the most damaging crater, the diameters and 
velocities subsequently calculated would represent upper bounds on ricochet debris sizes and 
speeds. Measured crater depths and diameters, as well as calculated crater volumes, for each of 
the 139 ricochet witness plates considered herein are presented in Appendix C. 

In the last phase of the analysis, equations for penetration depth and crater diameter were 
solved for particle diameter and velocity in terms of all other parameters, such as density, yield 
strength, wave speed, and so forth. Substitution of the appropriate parameter values in these 
equations yielded an estimate for the size and speed of the particle that produced a particular 
crater. This procedure was applied to the most damaging crater dimensions as identified 
previously. The penetration depth and crater mouth diameter equations are listed in Appendix D, 
some rewritten for consistency. The material property values used in these equations is presented 
in Table 5.1 below. 


Table 5.1 Material Property Values 


Symbol 

Property 

Value 

Units 

c„ 

Bumper Speed of Sound 

5.04 

km/s 

Pp 

Projectile Density 

2.718 


Pb 

Bumper Density 

2.718 

BjiMl 

Hb 

Brinell Hardness Number 

130 

KKSfl 

s b 

Bumper Dynamic Hardness 

6.37E+10 


S 

Bumper Shear Strength 

2.83E+09 

CEE51 

Svb 

■ A fvuminljEB f» E85BSM 

1.85E+10 


Y b 


2.78E+09 

E52EM 

B b 

Bumper Hardness 

1.27E+10 


Eb 

Bumper Elastic Modulus 

73.8 

GPa 


Since there are 12 penetration depth equations and 6 crater diameter equations, this 
method should have resulted in 72 estimates for the diameter and 72 estimates for the velocity of 
each crater-producing projectile. However, equations (D.l 1) and (D.12) were not used in 
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subsequent analyses because the upper limit of the velocity regime for which they are valid is 
much lower than that of the other penetration depth equations. Additionally, in the process of 
pairing the penetration depth and crater diameter equations, it became evident that not all 
equation pairs were compatible. Because of the exponential form of the equations, certain 
combinations of equations led to powers of zero for an unknown diameter or velocity. These 
particular equation pairs, therefore, could not be used to solve for the unknown quantities. This 
situation is analogous to finding the intersection of two parallel lines in Euclidean georretry. 
Specifically, penetration depth equations with a V 20 term could not be paired with crater diame ter 
equations having a V 2 term. Thus, in order to obtain unique solutions for particle velocity and 
diameter, depth equations with a V M term could only have been paired against diameter equations 
without a V 2 term, while depth equations without a term were paired against Hiarrv»tpr 
equations with a V 2 term. 

Furthermore, even though an equation pair did produce a solution, the resultant particle 
size occasionally exceeded that of the crater diameter, sometimes by a factor of three or four. 
However, it was previously shown that the heated material surrounding a high-speed inpact 
crater relaxes as it cools after the impact event, which can cause a reduction in crater diameter 
and depth of approximately 20-25%. Therefore, while it is possible that a crater could have been 
produced by a particle whose diameter exceeded the size of the crater opening, it is unlikely that 
the diameter of the particle could have exceeded the surface diameter of the crater it produced by 
more than 25%. As a result, a particle diameter value greater than 1.25 times a corresponding 
measured crater surface diameter was rejected. 

Figures 5. 1 and 5.2 show plots of equations (D. 1-D. 10), the penetration depth equations, 
and (D. 1 3-D. 1 8), the crater mouth diameter equations, as a function of inpact velocity for the 
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Normalized Penetration Depth (p/d) 


material parameter values given in Table 5.1. Examination of these plots reveals several 
interesting characteristics of the crater depth and diameter equations. 



-♦-Eqn(DI) 
-•-Eqn (D2) 
Eqn (D3) 
-*-Eqn(D4) 
-»-Eqn(D5) 
-*-Eqn(D6) 
—♦—Eqn (D7) 
-•-Eqn (D8) 
-•-Eqn (D9) 
-»- Eqn (DIO) 



-•—Eqn (D13) 
-•-Eqn (D14) 
-«-Eqn(D15) 
-•-Eqn (D16) 
-t*-Eqn(D17) 
-*-Eqn(D18) 
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consistent in their prediction trends actually values, however, cm vaiy 
significantly. Because of a lack of corroborative information for the trends and values 
predicted by equation (D.8), it was not considered in any of the subsequent analyses. 
2) The oato* equations appear to fall into two fairly distinct groups with regard to both 


predictive trends as well as predicted vetoes. Within each group, however, the 
predk&diratoes are fairly consistent. 


Based on these observations and the comments made previously regarding the pairing 


combinations were used to calculate candidate ricochet particle velocity-diameter values: 


Table S.2 Penetration Deg&-€!fatar Diameter Equation Pairs 



Velocity 

Twin 


Velocity 

Term 

D.l 

v *V 

D.14 

K WU 

D.2 

yM 

D.I4 

■BUN 

D.3 

yW 

D.14 

mmtk 

D.4 

yH 

D.14 

■sa 

D.5 

y®3 

D.14 

u 

D.7 

yW 

D.14 


D.10 

y»3 

D.14 

mam 

D.6 

Ujl | 

D.13 

; 

D.9 


D.13 

v M 1 


These considerations reduced the number of calculated ricochet particle velocity and 
diameter value pairs for each most damaging crater from 72 to 9 or less. The resulting 
particle diameters and velocities corresponding to the depths and diameters of the most damaging 
craters (taken from Appendix D) are given in App e nd ix E. In Appendix E, ‘Vx-y’ and ‘dx-y’ refer 
to the a particle velocity or diameter, respectively, calculated using a combination of crater depth 




















equation (D.x) and crater mouth diameter (D.y) from Appendix D. Grayed-out areas are 
calculated particle velocity-diameter combinations that are not valid, most likely because the 
calculated ricochet particle diameter exceeded the crater mouth diameter (indicated by a value of 
‘dx-y/d’ that is greater than one. 

For each test, valid particle velocity-diameter were reviewed to determine two max-min 
combinations for subsequent regression analyses: Vmax and the corresponding dnun. and V min and 
the corresponding dmax- In this manner, upper and lower bounds on velocity and size can be 
formed for the most damaging ricochet debris particle to be created in a given impact scenario. 
These max-min values are provided in Appendix F. 

Four empirical predictor equations for were developed using the data in Appendix F. 
These equations can be used to calculate V™*, 6 **, V™,, and <!«» in terms of bumper thickness 
and impact parameters, and were all in the following form: 


V: 


I 


= A 


(v V 

P 

M 

l c »J 

KJ 


cos D 0 p + E , i = max, min 


(5.1a,b) 


ir* 


f y N 8 '' \ c 


\ Cb J 


v dp 


cos D 0„ + E , i = max, min 


(5.2a,b) 


Table 5.3 below presents the values of the regression coefficients A-E and the correlation 


coefficients for equations (5.1a,b) and (5.2a, b). 

Table 5.3 Parameter Values and Correlation Coefficients for Equations (5. 1) and (5.2) 


Equation 

Quantity 

A 

B 

C 

D 

E 

Correlation 
Coefficient (R 2 ) 

5.1a 

Vmax 

0.4294 

-1.8335 

-0.2799 

-0.2562 

0.3384 

0.417 

5.1b 

dmin 

-0.6799 

-0.08769 

0.01119 

1.0558 

0.5998 

0.712 

5.2a 

Vmin 

0.3339 

-1.2209 

-0.1002 

-0.1588 

0.2206 

0.254 

5.2b 

dmax 

0.5732 

-0.02872 1 

-0.04935 

-0.4569 

-0.4978 

0.747 
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As can be seen from Table 5.1, the equations for dm and d** have reasonable R 2 values 
while those for dan and d** are somewhat low. This indicates that there is a fair degree of scatter 
in the ca lc u la ted ricochet particle velocity values, while the level of consistency in the calculated 
diameter values is Surly high. It is not clear at this time why this has occurred, especially since 
both velocity and diameter quantities were calculated simultaneously using the same data and the 
same equations. 

Equation (5. la,b) and (5.2a,b) can be used to obtain a bound on the velocity and diameter 
of the most dsrnagjiv ricochet debris particle that would be treated in a given oblique 
hypervelocity impact event. However, these equations must be paired appropriately: V*, must be 
paired with d^ while Vm must be paired witii (W This will provide, for example, upper and 
lower limits of expected ricochet particle velocity and the particle diameters corresponding to 
those velocities. 

Figures 5.3 through 5.6 below show plots of equations (5.1a,b) and (5.2a,b) for an initial 
projectile diameter of 0.795 cm, a 0. 127 bumper thickness, for impact velocities ranging between 
3 and 8 km/s, and for initial trajectory obliquities of 30°, 45°, 60°, and 75°. In these plots, the open 
tick marks represent values calculated using equations (51a,b) and (52a,b), while the solid tick 
marks represent simple numerical averages of corresponding calculated values. 

Table 5.4 below presents a comparison between the average ricochet debris particle 
diameters and velocities presented in Reference [2] and the average particle velocities and 
diameters calculated using equations (5. la,b) and (5.2a,b) under the same impact conditions. As 
can be seen from this table, the average diameter values predicted by the equations developed in 
this study compare favorably with those obtained previously. However, the average velocity 
values calculated using equations (5. la,b) and (5.2a,b) are approximately twice the values 
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reported previously. These differences and similarities serve to 1) reinforce the need to explore 
further the particle velocities obtained using the technique developed in this study, and 2) increase 
the confidence in the particle diameter values obtained using this technique. 

Table 5.4 Comparison of Average Ricochet Particle Diameters and Velocities 



d«* (cm) 


e. (deg) 

Reference [2] 

This Study 

Reference [2] 

This Study 

45 

0.174 

0.121 

2.07 

4.25 

60 

0.221 

0.204 

2.01 

4.42 

75 

0.357 

0.348 

2.35 

4.78 


d«r« (cm) 

(km/s) 1 

d„ (cm) 

Reference 121 

This Study 

Reference [21 

This Study 

0.475 

0.203 

0.164 

2.17 

4.35 

0.635 

0.258 1 

0.224 ! 

2.15 

4.49 

0.795 

0.303 

0.285 

2.08 

4.61 



Impact Velocity (tanfe) 

Figure 5.3 Va*, Vm»x and dam, dam for a 30° Impact 
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dmin, dmax (cm) 




VMn, VtMx (hmA>) 


•GhVtat 

-B-VMn 


0.00 1 I I ' I t I 1 i < i i 1 1 i I 0.000 

3.00 3.40 3.80 4.30 4.60 9.00 3.40 5 J 0 6.20 6.60 7.00 7.40 7.80 

jfinift Vainnjv ftnM 

Figure 5.4 V*„ and d«k, dUn for a 45* Impact 


-a-Vtmx .. 0.-I00 

-»-VWn 

-B-Vavg 

-O-dnwc . . 0.060 


1 1 

o.oo \ ■■■» \ 1 1 — — 4— < 1 — * — ♦ -4 - i > — L o.ooo 

300 3.40 3.80 4.20 4.00 9.00 9.40 5.80 8.20 6.80 7.00 7.40 7.80 

torpmdt VatooHy (taTVt) 

Figure 5.5 V.*, V«,* ami dU* <U« for a 60“ Impact 
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(u») nup ‘upup 





(•/UB0XMUA ’U)UIA 


10.00 T 


T aaoo 



+ 0.900 


4-0.400 



q 9 O > O o c 


9 9 9 OOC» 9 ft0 ^ ' 


1 0.300 



0.00 


-O-Vmw 

-B-Vrrln 

-*-V»vg 

-O-drmx 


+ 0.200 


+ 0.100 


I I , +- 1 1 1 1- 1 1 1 • — L 0.000 

3.00 3.40 3.80 4.20 4.80 5.00 5.40 5.80 6.20 6.80 7.00 7.40 7.80 

Impact Velocity (km/e) 


Figure 5.6 V™,, V« and <U, <L« for a 75* Impact 
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dmin, dmax (cm) 




€.0 SUMMARY AND RECOMMENDATIONS 

6.1 Summary 

An empirical model that that characterizes the secondary qjecia created by a high speed 
impact on a typical aerospace structural surface has been successfully developed. This model 
developed provides the following information as a function of impact parameters (speed, angle, 
projectile diameter) and target plate geometry (e.g. thickness, etc): 



cloud center-of-mass; 

• average velocity of the ricochet debris cloud material; and, 

• velocity and mass of the largest particle(s) in the ricochet debris cloud. 


The angles defining the spread of the ricochet debris cloud and the trajectory of the debris 
cloud center-of-mass were obtained using die spatial distributions of ricochet debris particle 
impacts on ricochet witness plates from over 200 high speed impact tests. The average velocity of 
the ricochet debris cloud is obtained using a model that characterizes the masses, trajectories, and 
velocities of the debris clouds created in an oblique high-speed impact. This model employs the 
three conservation principles, elementary shock physics theory, and fundamental thermodynamic 
principles to obtain a system of algebraic equations for the various debris cloud masses, 
trajectories, and velocities. Finally, relationships for crater diameter and depth are applied to the 
deepest craters in each ricochet witness plate to "back out" the diameters, masses, and velocities 
of the ricochet debris cloud particles that created these craters. This information is then used to 
develop empirical relationships that predict the velocity and mass of the largest ricochet debris 
cloud particle in terms of impact parameters and bumper plate thickness. 
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6.2 Recommendations 

Based on the work performed, the following recommendations are made for continued 
activities in this area. 

6.2.1 Ricochet Debris Qoud Spread Angie Modelling 

1) The discrepancy between empirical observations ami SPH predictions of ricochet debris cloud 
spr ea d should be explored and reconciled. It is suggested that an alternative means of defining 
ricochet debris cloud spread needs to be developed, one that will allow the successful use of 
empirical as well as numerical data. 

2) Thus far, ricochet debris cloud spread angle modelling efforts have focussed on characterizing 
the spread of the debris cloud particles “in the plane of the impact trajectory”. Future efforts 
should focus on the spread of the debris cloud out of this plane. 

6.2.2 Ricochet Debris Cloud Velocity Modelling 

Efforts should continue to reduce the dependence of the model on empirical or user- 
controlled parameters. A preliminary effort involving oblique shock wave theory was successfully 
completed by the author [8]; however, the modelling effort was at a level of complexity that is 
inconsistent with that employed in the model presented in this report. Some aspects of oblique 
shock wave theory should, however, be explored and considered for implementation in the debris 
cloud model presented herein. 

6.2.3 Ricochet Debris Qoud Particle Diameter and Velocity Modelling 

The reasons for the differences between the values of the correlation coefficients for the 
empirical diameter and velocity equations should be explored and reconciled, including the 
consideration of an alternative equation form should be considered for the ricochet debris cloud 
particle velocity. 
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APPENDIX A 


EMPIRICAL TEST PARAMETERS AND RESULTS 
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Table A-l Empirical Test Parameters and Results, Phase B NASA/MSFC Test Series 


Test 

tb 

dp 

v. 

Op 

0r 

£ . . V 

No. 

WB5M 

ITCIiDiH 

WLGDU 

■Gl» 

■rrriM 

m\±Jm 

001A 

0.203 

0.795 

6.62 

45 

11.3 

26.3 

OOIB 

0.203 

0.795 

6.53 

45 

12.6 

26.9 

002A 

0.160 

0.795 

6.50 

45 

9.5 

23.0 

002B 

0.160 1 

0.795 

6.45 

45 

12J_ 

26.0 

003 A 

0.102 

0.795 

6.54 

45 


20.9 

004A 

0.203 ! 

0.795 

6.28 

65 

mm 

25.3 1 

136A 

0.160 

0.635 

6.25 

55 

8.6 

25.1 

136B 

0.160 

0.635 

7.24 

55 

10.8 

29.2 

136C 

0.160 

0.635 

6.67 

55 

13.0 

28.8 

137D 

0.081 

0.635 

7.03 

45 

9.6 

24.7 

150A 

0.160 

0.635 

7.00 

45 

11.4 

24.2 

151A 

0.203 

0.635 

6.88 

45 

13.5 

21.5 

154 A 

0.102 

0.475 

6.83 

45 

11.1 

19.2 

155A 

0.160 

0.475 

7.02 

45 

12.6 

18.1 

156 A 

0.160 

0.475 

7.10 

65 

6.8 

16.9 

156B 

0.160 

0.475 

5.95 

65 

wtmam 

19.2 

156C 

0.160 1 

0.475 

4.15 

65 

7.1 

16.5 

157 A 

0.160 

0.475 

7.40 

60 

10.0 

22.3 

162 A 

0.160 

0.475 ' 

6.53 1 

30 

16.7 

31.0 

168 A 

0.081 

0.635 

5 54 

45 

8.8 

19.3 

168B 

0.081 

0.635 

5.98 

45 

10.8 

25.1 

168C 

0.081 

0.635 

6.67 

45 

19.4 

29.9 

168D 

0.081 

0.635 

7.02 

45 

21.7 

30.8 

169B 

0.081 

0.635 

6.55 

45 

10.8 

25.1 

201A 

0.102 

0.635 

4.33 

45 

10.4 

23.3 

201B 

0.102 

0.635 

5.51 

45 

8.8 

23.1 

20 ID 

0.102 

| 0.635 

7.59 

45 

14.3 

22.8 

202C 

0.102 

! 0.475 

5.25 

45 

10.8 

23.3 

202D 

0.102 

0.475 

6.44 

45 

6.2 

14.6 

202E 

0.102 

0.475 

7.19 

45 


16.0 

203 A 

0.102 

0.762 

4.79 

65 

■EZ3 

21.6 

203B 

0.102 

0.762 

3.65 

65 

■Q 

19.8 

203C 

0.102 

0.762 

2.72 

65 

! 6.8 

21.0 

203D 

0.102 

0.762 

5.59 

65 


22.2 

203E 

| 0.102 

0.762 

6.72 

65 

8.1 

23.5 

203F 

0.102 

0.889 

3.05 

65 

85 

26.0 

203G 

0.102 

0.889 

4.64 

65 

8.1 

25.5 

204A 

0.102 

0.635 

4.77 

65 

8.5 

24.2 

204B 

[ 0.102 

0.635 

5.86 

65 

mm 

24.0 

204C 

0.102 

0.635 

4.25 

65 

1 9.0 

28.2 
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227 B 

0.081 

0.635 

7.19 

45 

230 A 

QEna 

0.475 

4.41 

45 

230 B | 

0.160 

0.475 

3.23 

45 

23 0 C 

0.160 

0.635 

5.18 

45 

230 D 

0.160 

0.635 

5.55 

45 

230 E 

0.160 

0.635 

6.57 

45 

231 A 

0.160 

0.475 

3.34 

65 

23 IB 

0.160 

0.475 

2.44 

65 

23 1 C 

0.160 

0.795 

6.59 

65 

23 ID 

0.160 

0.795 

7.26 

65 

301 - 

0.160 

0.635 


45 

303 - 

0.160 

0.795 

4.65 

45 

303 A 

0.160 

0.795 

3.72 

45 

303 B 

0.160 

0.795 

4.42 

45 

306 - 

0.160 

0.953 

6.35 

45 

319 - 

0.102 

0.795 

2.99 

45 

320 - 

0.160 

0.795 

308 

45 

321 - 

0.203 

0.795 

3.01 

45 

324 - 

0.102 

0.795 

4.12 

45 

325 - 

0.160 

0.795 

msm 

45 

326 - 

0.203 

0.795 

1 4.25 

45 

333 - 

0.102 

0.475 


45 

334 - 

0.102 

0.475 

1 3.66 

45 

335 - 

0.102 

0.635 

msm 

45 

336 - 

0.102 

0 635 

4.54 

45 

336 A 

0.102 

0.635 

5.76 

45 

337 - 

0.102 

0.795 

6.90 

45 

338 - 

0.102 

0.795 

7.02 

45 

339 - 

0.102 

| 0.953 

IBlSi 

45 


13.8 

18.8 
10.8 
11.2 
10.4 
12.7 

5.4 

5.9 

8.1 




8.2 

8.5 

10.1 

8.4 

9.8 

8.4 

9.3 

8.7 

11.0 

9.6 

10.6 


9.9 

12.5 

12.2 

124 

8.5 


25.6 

23.9 

21.8 


27.0 

24.9 

19.4 

27.2 

20.9 

23.3 

19.3 

21.0 

17.5 

18.3 

18.4 

20.0 

22.6 
23.3 

24.0 
23.3 
25.2 

25.5 

24.0 

23.1 

20.9 

24.2 

25.0 

25.2 

26.9 

























































































































































Table A-2 Empirical Test Parameters and Results, Phase C/D NASA/MSFC Test Series 


4001-A 


4001-B 


4001-C 


4001-D 


0.203 


4002-B I 0.203 


4002-D 


4002-E 


4003-B 


4003-C 


4003-D 


4004-A 


B 


4004-C 


4100-A 


4100-B 


4101-B 


4101-C 


4101-D 


4102-A 


4102-C 


4102-C1 


4102-C2 


4102-D 


4103-A 


4103-B 


4103-C 


4103-D 


4104-A 


4104-B 


4104-C 


4104-D 


4105-A 


4105-A1 


0.127 


0.127 


0.127 


0.127 


0.127 


0.127 


0.127 


0.127 


0.127 


0.127 


0.127 


0.127 


0.127 


0.635 


0.795 


0.795 


0.795 


0.475 


0.475 


0.475 


0.635 


0.635 


0.795 


0.795 


6.14 

45 

9.0 

28.4 


12.8 

28.4 

mm 

45 

8.3 

TT.l 

4.12 

45 

5.2 

14.3 

6.24 

45 

■oa 

19.2 

6.05 

45 

11.3 

29.5 

6.02 

45 

6.3 

19.0 

7.18 

45 

6.4 

19.3 

2.94 

60 

10.4 

25.8 

3.98 

60 

mm 

25.1 

5.88 

60 

20.0 

37.4 

7.37 

60 

5.5 

13.7 

7.23 

60 

9.5 

29.1 

4.19 

60 

8.1 

27.6 

6.12 

60 

8.0 

28.9 

7.52 

60 

15.9 

29.5 

2.92 

60 

4.9 

21.0 

2.98 

60 

■o 

23.5 
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4105-B 

0.127 

0.795 

4.02 

60 

4105-C 

0.127 

0.795 

6.15 

60 

4105-D 

0.127 

0.795 

7.23 

60 

4106-A 

0.127 

0.475 

3.05 

60 

4106-A1 

0.127 

0.475 

3.10 

75 

4106-B 

0.127 

0.475 

4.12 

60 

4106-B1 

0.127 

0.475 

3.99 

75 

4106-C 

0.127 

0.475 

5.95 

75 

4106-D 

0.127 

0.475 

7.56 

75 

4107-A 

0.127 

0.635 

3.05 

75 

4107-B 

0.127 

0.635 

4.11 

75 

4107-C 

0.127 

0.635 

6.20 

75 

4107-D 

0.127 

0.635 

7.64 

75 

4 108- A 

0.127 

0.795 

3.12 

75 

4108-A1 

0.127 

0.795 

2.95 

75 

4108-B 

0.127 

0.795 

3.97 

75 

4108-C 

0.127 

0.795 

5.96 

75 

4108-D 

0.127 

0.795 

7.07 

75 

4109-A 

0.203 

0.475 

3.27 

45 

4109-B 

0.203 

0.475 ] 

4.14 

45 

4109-C 

0.203 

0.475 

6.53 

45 

4109-D 

0.203 

0.475 

7.46 

45 

4110-A 

0.203 

0.635 

3.25 

45 

4110-B 

0.203 

0.635 

4.00 

45 

4110-C 

0.203 

0.635 

5.76 

45 

4110-D 

0.203 

0.635 

6.96 

45 

4111-A 

0.203 

0.795 

2.85 

45 

4111-B 

0.203 

0.795 

3.94 

45 

4111-C 

0.203 

0.795 

5.97 

45 

4111-D 

0.203 

0.795 

6.81 

45 

4112-A 

0.203 

0.475 

3.33 

60 

4112-B 

0.203 

0.475 

4.05”| 

60 

4112-C 

0.203 

0.475 

5.87 

60 

41 12-D 

0.203 

0.475 

7.50 

60 

4113-A 

0.203 

0.635 

2.97 

60 

4113-B. 

0.203 

0.635 

! 3.77 

60 

4113-C 

0.203 

0.635 

I 6.30 

60 

4113-D 

0.203 

0.635 

7.12 

60 

4114-A 

0.203 

0.795 

3.13 n 

60 

4114-B 

0.203 

0.795 

3.98 

60 

4114-C 

0.203 

0.795 

5.92 

60 

4114-D 

0.203 

0.795 

7.40 

60 

4115-A 

0.203 

i 0.475 

3.13 

75 

4115-B 

0.203 

0.475 

4.08 

75 


7.3 

25.1 

m 

24.8 

m 

25.5 

9 4 

26.8 

4.8 

16.3 

12.7 

287 

4.1 

16.1 

■o 

19.0 

5.8 

22.4 

6.5 

KOI 

5.8 

24.9 

m 

20.6 

6.1 

23.0 

■eo 

22.6 

5.4 

18.7 

WKSM 

22.8 

6.9 

23.5 

6.5 

22.6 

14.7 

39.4 

15.3 

31.8 

13.2 

29.5 

25.6 

34.8 

16.1 

32.3 

22.3 

32.3 

18.4 

30.6 

17.5 

33.6 

19.3 

33.0 

17.6 

35.0 

20.9 

35.5 

18.4 

36.9 

10.2 

29.1 

8.9 

26.8 

12.8 

31.0 

14.4 

28.0 

7.1 

29.8 

12.1 

41.0 

10.1 

28.0 

10.2 

31.0 

9.5 

27.8 

12.4 

32.0 

13.3 

28.9 

10.0 

31.8 

■Bl 

19.5 

6.1 

13.1 


I- 
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4115-C 


| 0.475 

1 6.06 

75 

l 

30.0 

4115-D 

El WM 

WX5M 

7.30 

75 

HI 

24.7 

4116-A 

EMM 

■«JW 

2.92 

75 

4.4 

17.5 

4116-B 

EEHS 


4.48 

75 

8.7 

27.4 

mmsm 


fabM 


75 

58 

20.3 

4116-D 

EMM 

ET5M 

7.36 

75 

Kl 

■EEQ 

rnsnsm 

EEM 

EMM 

3.11 

■a 

5.3 

18.0 | 

mamm 

EMM 

0.795 

El 

75 

■EE1 

■Egg 

4117-C 

EMM 

EMM 

6.03 

75 

6.5 

26.9 

4117-D 

EMM 

0.795 | 

7.20 

75 

6.8 

25.5 
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Table A-3 Empirical Test Parameters and Results, NASA/MSFC EH Test Series 



EHRP3 


EHRP4 


EHRP5 


EHRP6 


EHRP7 


EHRP8 


EHRP9 


EHSS4C 



0.160 


0.160 


0.160 


0.160 


0.160 


0.160 


0.160 





Or 

Qff 

* . V 


mmm 

■gzjb 

■BUM 

MUjJm 

0.795 

6.93 

75 


31.1 

0.795 

6.91 

75 

■Q 

27.2 

WfSjiEM 

6.82 

75 

9.3 

29.8 

0.795 

7.01 

45 

15.5 

29.3 

0.635 

7.22 

75 

6.6 

27.5 

0.795 

7.17 

60 

10.6 

27.9 1 

0.475 

7.52 

75 

8.2 

25.1 

0.795 

7.16 

75 

7.9 

29.0 

0.795 

6.93 

60 

10.6 

24.2 

0.795 

6.85 

65 

8.7 

21.2 

0.795 

6.83 

45 

12.5 

27.7 

0.635 

7.71 

60 

10.8 

26.0 

0.635 

7.56 

65 

13.0 

30.1 

0.635 

7.63 

45 

11.5 

24.5 

0.475 

8.04 

60 

18.5 

30.8 

0.475 1 

7.39 

45 

11.9 

27.2 

0.475 

7.34 

65 

13.8 

34.1 

0.635 

5.58 

45 

9.9 

28.9 
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Table B-l SPH Numerical Simulation Parameters and Results 


Test 

t b 


No. 


mssmi 

SPH-01 

0.160 

0.635 

SPH-02 

0.160 

0.635 

SPH-03 

0.160 

0.635 

SPH-04 

0.160 

0.635 

SPH-05 

0.160 

0.635 

SPH-06 

0.160 

0.635 

SPH-07 

0.160 

0.635 

SPH-081 

0.160 

0.635 

SPH-09 

0.160 

wsmi 

SPH-10 

0.160 

0.953 

SPH-11 

0.160 

0.953 

SPH- 12 

0.160 

0.953 

SPH-13 | 

0.160 

0.953 

SPH-14 

0.160 

0.953 

SPH- 15 | 

0.160 

0.953 

SPH-16 

0.160 

0.953 

SPH-17 

0.160 

0.953 

SPH-18 

0.160 

0.953 

SPH-19 

0.160 

1.270 

SPH-20 

0.160 

1.270 

SPH-21 

0.160 

1.270 

SPH-22 

0.160 

1.270 

SPH-23 

0.160 j 

1.270 

SPH-24 

0.160 

1.270 

SPH-25 

0.160 

1.270 

SPH-26 

0.160 

1.270 

SPH-27 

0.160 

1.270 

SPH-28 

1 0.127 

0.795 

SPH-29 

0.127 

0.795 

SPH-30 

0.127 

0.795 

SPH-31 

0.127 

0.795 

SPH-32 

0.127 

0.795 

SPH-33 

0.127 

0.795 

SPH-34 

0.127 

1.113 

SPH-35 

0.127 

1.113 

SPH-36 

0.127 

1.113 

SPH-37 

0.127 

1.113 

SPH-38 

0.127 

1.113 

SPH-39 

0.127 

1.113 

J- — 
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MEASURED CRATER DEPTHS AND DIAMETERS 
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Table C-l Measured Crater Depths and Diameters, Calculated Crater Volumes 
Phase B NASA/MSFC Test Series 


Test 

No. 


OOIB 


002B 


201A 


20SA 


205C 


20SD 


205E 


206E 


206 F 


2UB 


211D 


212B 


2 16 A 


216C 


217A 


217B 


217C 


217D 


217E 


218A 


218C 


221B 


221 C 


226A 


226B 


227A 


227B 


230B 


230C 


230D 


230 E 


301- 


303- 


303A 


303B 


319- 


321- 


324- 


325- 


326- 


333- 


V, 

(km/s) 


6.56 


6.51 


4.33 


4.20 


5.30 


6.42 


3.15 


3.24 


6.42 


5.88 



0.762 


0.889 


0.795 


0.795 


0.795 


0.635 


0.635 


0.635 


0.889 


0.889 


0.475 


0.475 


0.635 


0.635 


0.635 


0.635 


0.475 


0.635 


0.635 


0.635 


0.635 


0.795 


0.795 


0.795 


0.795 


0.795 


0.795 


0.795 


0.795 


0.475 


0.042 

0.095 

0.033 

0.113 

0.026 

0.055 

0.037 

0.084 

0.018 

0.071 

0.018 

0.074 

0.023 

0.089 

0.012 

0.070 

0.037 

0.061 

0.058 

0.105 

0.052 

0.105 

0.031 

0.096 

0.058 

0.101 

0.045 

0.119 

0.036 

0.098 


Vol, p, 
(xlO* cm 3 ) (cm) 


0.020 0.043 


Vob p, 
(xlO ’cm 3 ) (cm) 




d, Vol, 
(cm) (xlO' 3 cm 3 ) 


0.025 


0.10 0.036 


0.10 0.029 


0.018 


0.040 


0.10 0.063 


0.10 0.020 


0.10 0.012 


0.08 0.004 


0.08 0.015 


0.08 0.102 


0.027 


0.036 


.024 


0.16 0.023 


0.16 0.043 


0.16 0.026 


0.16 0.062 


0.16 0.016 


0.16 0.040 


0.10 0.042 


0.20 0.071 


0.10 0.016 


0.16 0.026 


0.20 0.105 


0.10 0.020 
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Table C-2 Measured Crater Depths and Diameters, Calculated Crater Volumes 
Phase C/D NASA/MSFC Test Series 


Test 

No. 


4001A 


4001B 


400 1C 


4001D 


4002B 


4002D 


4002E 


4003A 


4003B 


4003C 


4003D 


4004A 


4004B 


4101A 


4101B 


4101C 


4102A 


4102C 


4102C1 


4103A 


4103B 


4103C 


4103D 


4104A 


4104B 


4104C 


4104D 


4105A 


410SA1 


410SB 


4106A 


4106A1 


4106B 


4106B1 


4106C 


4106D 


4107A 


4107B 


4107C 


4107D 


4108A 


V P 

(km/s) 


3.15 


4.29 


6.12 


6.71 


3.97 


7.14 


6.41 


3.43 


6.29 


3.18 


6.22 


3.19 


6.08 


3.14 


4.13 


6.14 


2.95 


6.24 


6.05 


2.94 


3.98 


5.88 


7.37 


7.23 


4.19 


6.12 


7.52 


2.92 


2.98 


4.02 


3.05 


3.10 


4.12 


3.99 


5.95 


7.56 


3.05 


4.11 


6.20 


7.64 


3.12 



Pi 

d, 

(era) 

(c«) 


0.795 


45 0.795 


45 0.795 


45 0.795 


75 0.795 


0.795 


0.795 


0.795 


45 0.795 


45 0.795 


45 0.795 


75 0.795 


0.795 


0.635 


45 10.635 


45 0.635 


45 0.795 


45 0.795 


45 0.795 


0.475 


0.475 


0.475 


0.475 


0.635 


60 0.635 


60 0.635 


60 0.635 


.795 


.795 


0.795 


60 0.475 


75 0.475 


60 0.475 


0.475 


0.475 


0.475 


75 0.475 


75 0.635 


75 0.635 


75 0.635 


75 0.795 


IQ 

IQ 

IQ 

IQ 


K2QG 

Kir 



IQ 

IQ 


0.13 0.129 0.186 


0.13 0.160 0.241 


0.106 0.214 


0.13 0.161 0.263 


0.13 0.238 0.174 


0.13 0.254 0.164 


.214 0.406 



0.021 0.085 


0.053 0.103 


0.045 0.138 


0.029 0.104 


0.194 0.282 


.500 0.226 


.500 0.204 


.041 0.086 


.052 0.065 


0.019 0.091 


.058 0.079 


.211 0.202 


.291 0.151 


0.029 0.065 


0.029 0.085 


.041 0.124 


0.031 0.084 


.042 0.088 


0.035 0.131 


0.006 0.042 


.134 


0.041 0.086 


0.053 0.118 


0.115 


0.049 


0.079 0.114 


0.051 0.124 


0.071 0.111 


0.057 0.117 


0.139 0.117 


0.052 0.102 


0.128 


0.044 















































































































































































































































































































































































































































































































































































































































































































































Table C-3 Measured Crater Depths and Diameters, Calculated Crater Volumes 

NASA/MSFC EH Test Senes 


Test 

No. 

v, 

(km/s) 

0, 

(deg) 

H 

H 

Pi 

(cm) 


Vol, 

(xlO 1 cm*) 

1 P* 

(cm) 

d* 

(cm) 

Vdj 

(xlO 3 cm 3 ) 

P) 

(cm) 

d, 

(cm) 

Vol, 

(xlO* 3 cm 3 ) 

EHAB 

6.91 

75 

0.795 

0.16 

0.615 

0.734 

| 13.011 


0.686 

6.801 

0.483 

0.566 

6.076 


7.22 

75 

0.635 

0.16 

PTT31 

0.650 

8.213 

0.361 

0.602 

5.138 

0.310 

0.445 

2.411 

EHPC 

7.58 

75 

0.475 

0.16 

0.386 

0.599 

5.439 

0.318 

0.447 

2.495 

0.345 

0.422 

2.413 

EHRP1 

6.87 

60 

0.795 

0.16 

0.140 

0.254 

0.355 

0.094 

0.241 

0.286 

0.117 

0.244 

0.365 

mza 

B21 

65 

0.795 

ram 

0.371 

0.632 

5.819 

0.229 

0.445 

1.781 

0.211 

0.445 

2.188 E 

1 

6.78 

45 

0.795 

0.16 

0.165 

0.368 

1.170 

0.150 

0.320 


0.135 

0.343 

0.832 

ussa 

7.65 

mm 

0.635 

0.16 

0.152 

0.279 

0.465 

P1FTO 

6.371 

1.168 

0.157 

0.328 

0.884 

| EHRPS | 

7.51 

65 

0.635 

0.16 

0.305 

0.528 

3.339 

0.330 

0.546 

3.863 

0.203 

0.411 

1.795 


fESM 

45 

0.635 

0.16 

0.097 

0.201 


0.114 

0.267 

0.426 

0.084 

0.211 

0.1% 

EHRP7 

7.98 

60 

0.475 

0.16 

0.323 

0.488 

3.021 

0.254 

0.396 

1.564 

0.203 

0.465 

2.298 

EHRP8 

7.34 

45 

0.475 

0.16 

0.155 

0.262 

0.418 

0.137 

0.279 

0.558 

0.168 

0.295 

0.574 


7.29 

65 

0.475 

0.16 

0.108 

0.221 

0.276 

0.096 

0.266 

0.356 


0.197 

0.181 

|EHSS4C| 

5.53 

45 

0.635 

0.16 

0.078 

0.126 

0.049 

0.074 

0.139 

0.056 

0.045 

0.088 

0.014 
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APPENDIX D 


EMPIRICAL DEPTH AND DIAMETERS EQUATIONS 



Penetration Depth Equations 


p/d = 2.28(pp/pb AVp/Cbf 3 

[9] 

(D.l) 

p/d= 1.96(pp/pb) l/2 (V p /Cb) 2/3 

[10] 

(D.2) 

p/dp= 1.50(Pp/pb) 1/3 (PpVp 2 /2S b ) ,/3 

[11] 

(D.3) 

p/ d p = 2.35(pp/ P b) 0 ' 70 (Vp/Cb) 2/3 

[12] 

(D.4) 

p/dp = 0.63(ppVp 2 /Syt) ,/3 

[13] 

(D.5) 

p/dp - 0.428( P p/pb)° 537 (Vp/Cb)° 376 (Y b /pbCb 2 )'° 235 

[14] 

(D.6) 

p/dp = 8.355xl0‘ 4 pp2/3pb-l/3(V p 2 /Hb) 1/3 

[15] 

(D.7) 

p/dp = 2.00( Pp / Pb) 4-52 (Vp/Cb) 1/136 

[16] 

(D.8) 

p/dp = 0.3 1 1(pp/Pb)° ,7 (PpVp 2 /S b ) 1/3 

[17] 

(D.9) 

p/dp = 0.36(pp/p b ) 2/3 (PpVp 2 /Bb) ,/3 

[18] 

(D.10) 

p - 2.973x1 0' 7 dp U H b ' o25 pp 05 pb J0167 Vp 4/3 

[19] 

(Dll) 

p = 1.129x10%* 036 Hb-° 25 Pp° V 

[19] 

(D.12) 

Crater Mouth Diameter Equations 



ad 2 p/dp 3 = 34(pp/pb) 3/2 (Vp/Cb) 2 

[20] 

(D.13) 

ad 2 p/dp 3 = 0. 120(pp/pb) ,/2 (ppVp 2 /Sb) 0 '* 45 

[17] 

(D.14) 

ad 2 p/dp 3 = 30.25(pp/pb) 3/2 (Vp/C b ) 2 

[10] 

(D.15) 

ad 2 p/dp 3 = 44. 10(pp/pb) 3/2 (Vp/C b ) 2 

[12] 

(D.16) 

ad 2 p/dp 3 = 2.65xlO J9 p p 7A5 Pb' ,/J Vp 4/3 / Hb 

[15] 

(D.17) 

ad 2 p/dp 3 = O.l^pp/pb^PpVp 2 ^ 

[18] 

(D.18) 


where a=0.75 if p>d/2 and a=1.00 if p<d/2. 
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APPENDIX E 


CALCULATED RICOCHET PARTICLE VELOCITIES AND DIAMETERS 
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■ Jill 


M 


Test No 


EHAB 


V7-14 I D7-14 ID7-14/dH V10-14 I DIO-141 D10-14AI 


■CHSESHJ 
■CEniEEia 


EHPC 


EHRP1 0.117 0.244 


EHRP2 0.371 


■eei 

'ESZlCEEIEEn 




3.8281 0.237 


3.9491 0.0941 0.387 


2.560 0.178 0.483 


7.520 0.066 


0.1241 0.336 


EHRP6 0.114 


EHRP7 


EHRP8 10.137 


EHRP9 


■ Egg 

PlPTIPlPf! 


I =t \ t-i-'I 


001 B 0.117 0.257 


002B 


IEE23EE3 


EED MHMPESO 


0.152 

0.568 


0.100 

0.358 




0.180 


205D 0.074 0.234 


205 E 


0.178 


0.155 


21 IB 0.147 


211D 0.132 


i eeb I 


2.833 0.1181 0.459 


1.1191 0.1331 0.738 


3.2291 0.0771 0.429 


m 


8.766 ! 0.0691 


0.465 


5.3941 0.0821 


2.1301 0.093 




:v . ». : &J&} 




1.989 0.205 


2.531 0.266 




4.5491 0.176 


1.1161 0.2621 0.742 


ySvKWX 


• • . ■ 








217A 


217B 0.196 


iMj rgn 

|n»iyjnPT?i 

iji^n FTgp 

■■■EE3I 


4.863 

0.093 

0.348 




2.360 

0.168 

0.504! 

I«Ar4»I-1 

0.232 

0.933! 




3.649 

0.078 

! 0.403! 


g nroi ijrw 
lEESIEEISI 
■MEEST 

m usm 

p t'I's'ilFTm 

i ifeirrm 


21 7E 


218A 


8C 10.160 


0.051 
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APPENDIX F 


CALCULATED RICOCHET PARTICLE MAX— MIN COMBINATIONS 
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Table F-l Calculated Ricochet Particle Max-Min Combinations 
Phase B NASA/MSFC Test Series 
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Table F-2 Calculated Ricochet Particle Max-Min Combinations 
Phase C/D NASA/MSFC Test Series 
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Table F-3 Calculated Ricochet Particle Max-Min Combinations 
NASA/MSFC EH Test Senes 
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